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We used ab-initio quantum chemical methods, treating explicitly the strong correlation effects 
within the cobalt 3d shell, as well as the screening effects on the effective integrals, for accurately 
determining on-site and nearest-neighbor (NN) interactions in the Nao.3sCo02 — I.3H2O super- 
conducting compound. The effective ligand field splitting within the t^g orbitals was found to be 
8 ~ 300 meV, the a\ g orbital being destabilized compared to the e' g ones. The effective Hund's 
exchange and Coulomb repulsion were evaluated to Jh ~ 280 meV and U ~ 4.1-4.8 eV for the ai g 
orbitals. The NN hopping parameters were determined within the three t^g orbitals and found to be 
of the same order of magnitude as the t2 g ligand field splitting. This result supports the hypothesis 
that a three band model would be better suited than a one-band model for this system. Finally we 
evaluated the NN effective exchange integral to be antiferromagnetic and J = —66 meV. 



I. INTRODUCTION 

The layered cobalt oxides have been attracting a lot of 
attention for the last few years. This interest is driven 
by the remarkable properties of the Na x Co02 compounds 
and more specially their hydrated counterpart. Indeed, 
superconductivity was discovered in the Nao.3sCo02 — 
I.3H2O 1 , for the first time in layered oxides, beside its 
discovery in the cuprates. 

The host material is composed of C0O2 layers, where 
the cobalt atoms are located in distorted edge-sharing 
octahedra forming a two-dimensional triangular lattice. 
The sodium atoms are located in a plane in between 
the C0O2 layers. In the hydrated superconducting mate- 
rial, the water molecules are intercalated in between the 
sodium and the cobalt oxide layers. The water interca- 
lation is fragile and the system looses its water out of 
a hydrated atmosphere. It seems however, that at low 
temperature, the water and sodium cations order in a 
two-dimensional super-cell and adopt a frozen local ge- 
ometry similar to clusters of Na + ions embedded in ice2>. 
Despite this possible ordering, it is believed that the ef- 
fect of the water molecules is only sterioi. Indeed, the 
water inclusion induces a large separation of the C0O2 
layers, responsible for essentially uncoupled cobalt layers 
and a two-dimensional physics 3 . This two-dimensional 
character is assumed to be necessary for superconductiv- 
ity to occur. 

A simple formal charge analysis shows that the cobalt 
ions are 3.65+ for x = 0.35, that is about one third of 
Co 3+ ions and two third of Co 4+ ions. Wet-chemical 
redox analyses revealed however a cobalt oxidation num- 
ber somewhat lower— : 3.46+, that is closer to half Co 3+ , 
half Co 4+ ions. The distortion of the C0O6 octahedra 
observed in the superconducting material corresponds to 
a compression along the x + y + z axis of the Co coor- 
dination octahedron. This trigonal distortion induces a 
lowering of the Oh local point group to a c?3 subgroup, 
thus splitting the t^g orbitals in a a\ g and two e' g ones 
(see figure 0] in the results section) . Authors however 



disagree on the relative energies of the orbitals and the 
amplitude of the splitting. While some authors^ support 
the idea that the a\ g orbital is lowered compared to the 
two e' g , other authors come to the opposite conclusion^. 
Density functional Theory (DFT) calculations! agree on 
the fact that the a\ g band is less filled than the ones orig- 
inating from the e' g orbitals. However, the relative ener- 
gies of the different atomic configurations are not directly 
accessible to DFT calculations. Indeed, electronic corre- 
lation is assumed to be very strong in this system and 
can be expected to strongly influence the local atomic 
excitation energies between the different configurations 
associated with the hole localization either on the a\ g or 
on of the two e' g cobalt atomic orbitals. The question 
of the relative energy of these different configurations is 
however crucial in order to determine the pertinent de- 
grees of freedom to be taken into account in a simple 
model, able to describe the low energy properties of the 
system. Indeed, assuming that the a\ g orbital is much 
higher in energy that the e' g ones, one should naturally 
conclude that the pertinent model for the description of 
the superconductivity is a one band t — J type of model. 
Assuming that the a\ g orbital is now much lower in en- 
ergy that the e' g ones, one comes to a two-band model, 
while if the three orbitals are only weakly split, the per- 
tinent model should consider all of them at the time. No 
consensus is reached nowadays in the literature and there 
is a large controversy on the pertinent model to consider. 

The aim of the present paper is to determine the lo- 
cal orbital energies and effective coupling parameters 
between the cobalt 3c? orbitals. For this purpose we 
used embedded clusters calculations and quantum chem- 
ical ab-initio methods treating exactly the correlation ef- 
fects within the 3d shell as well as the screening effects 
that renormalize the interactions. Such methods allow 
the direct computation of the local parameters such as 
the atomic effective ligand field splitting, the one-site 
Hubbard U coulombic repulsion as well as the Hund's 
exchange. In addition, nearest neighbor exchange and 
transfer interactions can be directly and accurately com- 
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puted. 

The next section will shortly describe the method, sec- 
tion three will relate the results and finally the last sec- 
tion will be devoted to discussion and conclusion. 



II. METHOD AND COMPUTATIONAL 
DETAILS 

The method used in this work (CAS+DDCI&) is a 
configurations interaction method, that is an exact di- 
agonalisation method within a selected set of Slater 
determinants, on embedded crystal fragments. This 
method has been specifically designed to accurately treat 
strongly correlated systems, for which there is no single- 
determinant description. The main point is to treat ex- 
actly all correlation effects and exchange effects within a 
selected set of orbitals (here the 3<i shell of the cobalt 
atoms) as well as the excitations responsible for the 
screening effects on the exchange, repulsion, hopping, etc. 
integrals. 

The CAS+DDCI method has proved very efficient to 
compute, within experimental accuracy, the local interac- 
tions (orbital energies, atomic excitations, exchange and 
transfer integrals, coulomb repulsion etc.) of a large fam- 
ily of strongly correlated systems such as high T c cop- 
per oxides^, vanadium oxidesiS, nickel and cuprate fluo- 
rides^, spin chains and ladders^, etc. 

The clusters used in this work involve either one cobalt 
(CoOg) or two cobalt atoms (C02O10) and their oxygen 
first coordination shell (see figure . These fragments 
are embedded in a bath designed so that to reproduce on 
them the main effects of the rest of the crystal ; that is 
the Madelung potential and the exclusion effects of the 
electrons of the other atoms of the crystal on the clusters 
electrons. 

The electrostatic potential is reproduced by a set of 
point charges located at the atomic positions. The 
charges are renormalized next to the bath borders 
in order to obtain an exponential convergence of the 
Madelung potential. The convergence accuracy was set in 
the present work to the mili-electron-Volt. The method 
used for this purpose is a generalization^ of the Evjen's 
method—. The nominal atomic charges used in this work 
are the formal charges, that is +3.65 for the cobalt atoms, 
—2 for the oxygen atoms and +1 for the sodium atoms. 
The sodium atoms being located at two crystallographic 
sites, with fractional occupations, we renormalized the 
associated charges with the crystallographic occupation, 
thus using a mean field averaging of the Madelung po- 
tential. At this point we would like to shortly discuss 
the question of the Co valency and insertion of [H30] + 
ions. Indeed, it is clear that both the Co 3+ and Co 4+ 
atomic configurations of the Co atom are present in the 
system, however the actual average charge supported by 
the cobalt is under debate. A modification of the cobalt 
average charge would act on our calculations through a 
modification of the electrostatic potential seen by the 



computed cluster. A global shift of the electrostatic po- 
tential seen by the cluster would of course have no ef- 
fects on our results. However, a relative shift of the elec- 
trostatic potential seen by the cobalt and oxygen atoms 
would have a strong effect. Indeed, a reduction of the 
potential difference between the bridging oxygen orbitals 
and the cobalt magnetic orbitals can be expected to in- 
crease the effective exchange integrals through the lower- 
ing of the ligand-to-metal charge transfer configurations 
that mediate the interactions. On the contrary, an in- 
crease of the potential difference between these orbitals 
would decrease the effective exchange. 

The exclusion effects are treated using total ions 
pseudo-potentials 15 (TIP) on the first shell of atomic sites 
surrounding the clusters. 




FIG. 1: a) C0O6 and b) C02O10 clusters used in the present 
calculations. 

In this description the water molecules are not explic- 
itly considered. Indeed, we supposed correct the usual 
idea that their role is only the increase of the c parame- 
ter, uncoupling the C0O2 layers. This point is correctly 
treated in our calculations since the atomic positions 
(both in the quantum clusters and in the baths) are taken 
as given in the crystallographic data. One could argue 
that the water molecules should be taken into account 
in the clusters environment, however since water do not 
present a total net charge but only a dipole moment, its 
effect is expected to be very small. We thus neglected 
this aspect. 

The calculations presented in this work are specific for 
the hydrated compound. In the dehydrated Nao.3sCo02 
phase, the atomic positions differ from the one used here, 
and more specifically the C0O2 inter-layers distances are 
strongly reduced. Even if the atomic fractional positions 
were not changed by the dehydration, the electrostatic 
potential seen by the cobalt and oxygen atoms of the 
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C0O2 layers would be strongly modified and thus the 
different orbital energies (and other parameters) would 
be affected. Whether one would like to perform calcu- 
lations for other sodium concentrations, adequate struc- 
tural data as well as cobalt average valency should be 
taken. 

The calculations were done using the MOLCASi£ and 
CASDlii set of programs. The basis sets used can be 
found in reference^. The structural parameters were 
taken from the Nature paper of Takadai. 

III. RESULTS 

As mentioned in the previous section we performed two 
sets of calculations. The first one aimed at determin- 
ing the one-site effective parameters such as the cobalt 
3d energy splitting due to the ligand field. This orbital 
splitting is renormalized by the correlation effects within 
the 3d shell as well as by the screening effects due to 
the virtual excitations. The second type of calculations 
aimed at determining the Co-Co interactions, transfer 
between the different 3d, former t2 g , orbitals as well as 
the effective exchange integrals. 



A. The effective on-site 3d energy splitting 

Calculations with a formal Co 4+ cation were performed 
on the embedded C0O6 cluster. Figure El reports the 
first excitation energies and the dominant term of the 
associated wave functions. Figure [21 shows the Co atomic 
orbitals. 

One sees immediately that the a lg orbital is a 3d c 2_ r 2 
orbital in the a, b, c crystallographic axes. As stated in 
many papers it can be written in as 

d c 2_ r 2 — (d X y -\- d xz -\- dyz^) /"n/3 

if x, y, z are the Co-0 nearly orthogonal directions of the 
C0O6 octahedron. As expected the three low energy or- 
bitals are pointing between the oxygen atoms, while the 
high energy ones are directed toward the ligands. An- 
other important point to notice is the amount of Co(3d)- 
0(2p) hybridization in the cobalt orbitals. Indeed, while 
the two e' g orbitals do not present noticeable derealiza- 
tion on the neighboring oxygen ligands, the a\ g orbital 
exhibits some mixing with the oxygen 2p z orbitals - 
namely about 2/3 Co (3d) and 1/3 0(2p) — and the two 
e g orbitals are strongly hybridized with the oxygen 2p 
with for the two e g about 55% Co(3d) and 42% on the 
0(2p). 

Figure |21 shows us that in the ground state, the hole is 
located on the a\ g orbital, which is destabilized compared 
to the e' g . Figure 0] pictures this result as the effective lig- 
and field splitting associated with the trigonal distortion. 
This result is in agreement with the finding of the differ- 
ent DFT ab initio calculations. It however disagrees with 



the ligand field analysis of reference^. It is well known 
that when the transfer (overlap) between a 3d metal or- 
bital and the occupied ligand orbitals increases, the metal 
3d orbital is destabilized. This can be seen in a simple 
tight binding picture between a metal 3d and ligand 2p 
orbitals. At the second order of perturbation, the metal 
3d orbital is destabilized by the quantity t' 2 /(ed— £ P ) while 
the ligand 2p orbital is stabilized by the same value, due 
to their hybridization. This is presently the case, since 
the apex oxygens approach the cobalt plane in the dis- 
torted tetrahedron, and thus the overlap between the a\ g 
cobalt orbital and the oxygen 2p orbitals is slightly aug- 
mented. The computed resulting effective orbital energy 
splitting is of the order of 300 meV as can be seen on 
figure|31 that is in global agreement with — even if some- 
what larger than — the LDA+U estimation^ (0.2 eV) 
from the top of the a\ g and e' g bands. Finally, one should 
notice a further splitting between the e' g orbitals them- 
selves. This very small splitting is due to the fact that 
the electrostatic field generated by the sodium cations on 
the cobalt does not present a perfect three-fold symmetry. 
This value is however very small and can be neglected for 
any practical purpose. 

The a\ g -e g splitting can be extracted from the Co 
atomic d — > d excitations of higher energies. Typically 
the S — 3/2 and 5 = 5/2 states should be computed. 
We found a a\ g -e g splitting of 1.8eV and thus a e' g ~e g 
splitting of 2.1eV, in global agreement with the LDA eval- 
uations of 2.5eV— . 



B. The inter-atomic interactions 

As mentioned in the introduction, one of the open 
questions in the modeling of the present compound is 
whether only the a\ g band is important for the physical 
properties or whether one should consider a multi-band 
model. Indeed, Lee et proposed from LDA+U cal- 
culations that a crossover occurs between a single band 
behavior and a three band behavior as a function of the 
band filling. For an equivalent sodium concentration 
x < 0.5 they suggest a three band model and for x > 0.5 
a single band one. On the model point of view, while 
several authors argue, on analytical as well as numerical 
results, that proper superconductivity behavior cannot 
be found using a single band t — J type modeSi, other 
authors found superconducting pairing within a single 
band t — J model^. 

In order to address this question one should be able to 
accurately evaluate the effective transfer and exchange 
integrals between the three former t2 g — ► 2e' g + a\ g 
orbitals. It is crucial in their evaluation to properly 
take into account all Coulomb repulsions, exchange and 
screening effects, as well as the metal-ligands charge 
transfers. 

These effective transfer integrals between two nearest 
neighbor cobalt atoms can be extracted from the first 
electronic excitations of an embedded C02O10 cluster 
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FIG. 2: a) Co 4+ orbitals in the Nao.3sCo02 — I.3H2O compound, b) Label of the different atoms. The c axis is orthogonal to 
the figure plane. 

with Co 3+ - Co A+ mixed cobalt valency. The computed teraction between the cobalt a\ g and e' g 3d orbitals and 

lowest six states can be associated with symmetric and 

antisymmetric combinations of the atomic states shown 

in figure |3J As mentioned earlier, the effective transfer 

integrals are strongly mediated by the oxygen 2p orbitals. 

Figure shows the oxygen 2p orbitals bridging the Co - 

Co interactions. One can write the following matrix in- 
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FIG. 4: Schematic representation of the trigonal distortion 
ligand field splitting as found from the ab initio correlated 
calculations. Angles in the distorted structure correspond to 
the atomic positions given in reference-. 



FIG. 5: Oxygen 2p orbitals mediating the interactions be- 
tween the aig and e! g orbitals of the two Co atoms. One sees 
that the p z and p2 orbitals overlap with both the ai fl and 
e' g 2 cobalt orbitals, while the pi overlap with the e' gl cobalt 
orbitals. 
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where the I and r superscripts are associated with the 
two cobalt atoms, the £j diagonal energies are the effec- 
tive orbital energies, U Pj are the cobalt 3di-oxygen 2pj 
transfers and the tfj are the direct transfer integrals be- 
tween the ?>di and 3dj orbitals of the two cobalt atoms. 
The direct integrals are small, however non negligible due 



to the short Co-Co distance and the p-d hybridization 
as far as the a\ g orbitals are concerned. In fact the ^2 
direct transfer could be omitted, since this one is really 
very small. In addition, we will see later that we cannot 
explain our results without explicitly considering at least 
the £q 2 term. Such a 3c? — 2p model may be considered as 
a bit too complex for practical uses. In addition, a rapid 
analysis of the computed wave functions shows that the 
explicit contribution of the 0(2p) — > Co(3d) excitations 
are quite small, even if very important for the mediation 
of the interactions between the two cobalt atoms. In- 
deed, the weight of these configurations is less than 4% 
in the wave functions. We can thus reduce the previous 
matrix into an effective Hamiltonian on the sole cobalt 3d 
orbitals. All the effects of the oxygen 2p orbitals should 
however be taken into account to properly describe the 
physics. It results the following effective Hamiltonian, 
where both inter-atomic and intra-atomic coupling terms 
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appear between the a\ g and e' g2 orbitals. 



with Ai = ei - e p + 5(U - 2J H ) - U p 
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where Uj are the effective resulting transfer integrals (di- 
rect plus mediated by the oxygen ligands) between the 
3di orbital of one cobalt and 3dj of the other, tpo2 is the 
intra-atomic a\ g -e' g2 effective transfer resulting from the 
interactions with the oxygen 2p orbitals. This last tp$ 2 
integral is in fact quite surprising since one does not ex- 
pect such an intra-atomic effective transfer to take place. 
It however can easily be explained in perturbation theory. 
FigureElpictures the mechanism responsible for the effec- 
tive transfers between the a\ g and e' g2 orbitals of the two 
cobalt atoms, mediated by one oxygen orbital. The fact 
that both the a\ g and e' g2 orbitals present non negligible 
transfer integrals toward the same oxygen 2p orbitals re- 
sult in the fact that four different configurations couple, 
at the second order perturbation theory, to the ground 
state of the Co 4+ -Co 3+ ions. Indeed, one has the follow- 
ing expression for the effective transfers between the a± g 
and e' g2 orbitals of the Co atoms 



'on 



{top J 2 (t 0P2 Y 



^20 — 



"-00 



L 20 



A 

top? t 2Vz 



A 

^0p2^2p2 



hp z t 2 p z 



^0p2^2p2 



A 



02 



102 



(2) 
(3) 
(4) 



5(17 — 2Jh) — Up, the A 2 denom- 
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inator being in the Descloiseaux acception of the quasi- 
degenerate perturbation theory 2 ^ 
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with A 2 = e 2 - e p + 5(U - 2J H ) - U p . 
U is the Coulomb repulsion for two electrons in the same 
Co 3d orbital, Jh is the Hund exchange, and U p is the 
Coulomb repulsion in the oxygen orbitals. Finally, the 
last term corresponds to a renormalization of the a\ g or- 
bital energy 
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Similarly the perturbation theory yields, at the second 
order, for t%i and £22 
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If the tp 2 o term is not included in the effective model, 
it is impossible to fit the computed wave functions and 
energies with a good accuracy. Of course, it is possible to 
rotate the orbitals in such a way to nullify the tp 2 Q term. 
It would however result in a set of orbitals that would 
be i) bond dependent and ii) no more belonging to the 
irreducible representation of the whole system. Indeed, 
the axis of the a\ g orbital would for instance be tilted 
with respect to the c axis in a direction dependent of 
the Co-Co direction. In order to keep a simple model, 
respecting the system symmetry, it appears better to us 
to use a non-zero tp 2 o term. 

The quality of the model can be evaluated by two ways. 
The first criterion is the norm of projection of the com- 
puted wave functions on the configurations space associ- 
ated with the model (model space). In the present case 
the model space is the six combinations of the atomic 
electronic states pictured in figure Large norms war- 
rant that the space supporting the model captures the 
physics of the system. In our calculations the minimal 
value obtained for the norms of the six projected wave 
functions is 0.87. The missing part of the computed wave 
function is composed by the excitations responsible for 
the screening effects that is more than 17 millions config- 
urations in our calculation. Table [I] displays the effective 
transfer parameters obtained from the fit of Hd- The 
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-89 -53 



TABLE I 

t2g — * dl 



Effective hopping parameters between former 
+ 2e' g orbitals. 



present fit was done by fulfilling (at the best) the stan- 
dard effective Hamiltonians requirements 2 ^ ; that is 

• that the projection of the exact wave functions in 
the space supporting the effective Hamiltonian are 
eigenf unctions of the later, 

• associated with the exact energies. 

It results in the minimization over the Hd parameters of 
the following Lagrangian 

£2 = \ H dP^exact{m) - E exact (m) exact M | 2 



where P is the projection over the space supporting Hd 
and ^ exact (m) are the ab initio wave functions associated 
with the ab initio energies E eX act(jn). The quality of 
the fit is very good since the average error (£/6) can be 
evaluated to lmeV. 

Let us now concentrate on the effective exchange in- 
tegral J between the a\ g Fermi level orbitals. J can be 
evaluated from the singlet-triplet excitation energy on a 
Co 2 Oio embedded cluster with two formally Co 4+ ions. 
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FIG. 6: Through bridge perturbative mechanism responsible for transfers between the a\ g and e' g2 orbitals of two nearest 
neighbor cobalt. The representation of the left and right Co atoms as well as the orbital denominations are given in the starting 
step. The arrows picture at each step the electron transfer yielding the present configuration. Under each configuration of the 
last step, the hopping or energy term to which it contributes is indicated. Analytic formula associated with the present picture 
can be found in the set of equations [2] to 



Our calculations yield an antiferromagnetic coupling, in 
agreement with the experimental findings 25 and LDA+U 
calculations for large U values^. The computed value is 

J = — 66meV 

One should note that our evaluation of both the ef- 
fective transfer and exchange integrals are larger than 
the LDA and even LDA+U evaluations found in the lit- 
erature. This fact is due to the well known problem of 
the density functional methods to correctly treat strongly 
correlated systems. Indeed, it is well known that in nu- 
merous cases, the DFT results strongly overestimated 
the ferromagnetism, in particular in transition metal ox- 
ides where 3d orbitals play an important role. This is 
in particular the case for the present compound since 
LDA finds, for all compositions from x = 0.3 to x = 0.7, 
the system ferromagnetic and metallic, in disagreement 
with experimental results. LDA+U somewhat correct 
this problem, however only in a mean-field way, the quan- 
tum fluctuations due to the electronic correlation being 
ignored. In the present case, the authors of reference^ 
show that the LDA+U method of incorporating correla- 
tion effects is ill-suited for the Na x Co02 family of com- 
pounds. 

C. One site bi-electronic repulsion and Hund's 
exchange 

The effective Coulomb repulsion of two electrons in the 
same a\ g orbital, U, can be extracted from the present 
calculations in two different ways. On one hand, U 
can be evaluated from the ground-state energy differ- 
ence between the Co 3+ and Co 4+ ions, embedded in the 



Na .35CoO 2 - 1.3H 2 0. It yields U x = 4.1 eV. On the 
other hand, U can also be evaluated from the knowledge 
of both the hopping and exchange integrals between two 
nearest neighbor a\ g Co orbitals. Its value can thus be 
estimated to Uz = 4.8 eV. One immediately notices that 
the first estimation is somewhat weaker than the second 
one. Let us analyze this discrepancy. 

• The first value (U±) corresponds to a static increase 
of the cobalt charge, while the second (U2) corre- 
sponds to the energy of the quantum fluctuations : 
Co 3+ -Co 5+ in a globally Co 4 +-Co 4 + state. The 
Hubbard and related models use a unique param- 
eter for these two concepts. As far as the raw re- 
pulsion integrals are concerned, there is indeed a 
unique U. However, screening effects act differ- 
ently on the static Co 3+ configuration leading to 
U\ and on the quantum Co 3+ -Co 5+ fluctuations 
leading to U^- Indeed, it is well known in quantum 
chemistry that the static screening acting on U\ is 
quite larger that the dynamical one acting on /7 2 - 
It thus results in a smaller value for U\ compared 
to U 2 - 

• The U2 value corresponds in fact to the difference 
between the one site repulsion U and the repulsion 
V between the a\ g orbitals of nearest neighbor Co 
atoms. It means that the on-site repulsion respon- 
sible for the correct quantum fluctuations between 
two nearest neighbor is even larger than the com- 
puted U2 value. The 4.86^ value thus constitutes 
a lower bound for the repulsion responsible for the 
quantum fluctuations. 
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Even-though one cannot definitely conclude on the rel- 
ative values of U and V, one can expect from the 
above analysis that the effective repulsion between near- 
est neighbor Co sites is probably small. The above U 2 
values were extracted for a one-band model. Whether one 
would like to use a three band model based on both the 
aig and the e' g orbitals, one should take off the screening 
effects due to latters. This can be done in a perturba- 
tive manner. Indeed, the difference in the screening of 
the U 2 parameter between the one-band and three-bands 
models can be evaluated to 

U 2 (lb) = U 2 (3b) + A 

where U 2 (lb) is the one-band U 2 and U 2 {3b) is the three- 
band one and 

t 2 t 2 t 2 t 2 J 2 

A = -2 — + 2-^-2 20 +2 2-^ 

S-2J H 5 S-U 8-U -2J H 26 

The on-site repulsion between two different 3d orbitals 
Udd' was taken in the classical way equal to U — 2Jh- 
Using the numerical values of the present work we find 

A = -I93meV thus Z7 2 (36) = 5.0eV 

Finally we would like to recall that the Coulomb re- 
pulsion in the a\ g orbitals was evaluated from soft X-ray 
photo-emission spectroscopy 2 ^. The cobalt core 2p spec- 
trum exhibits well separated Co 4+ and Co 3+ levels, con- 
sistent with a repulsion value in the range of U ~ 3-5 eV, 
in total agreement with our calculations. Let us note that 
these values are somewhat weaker than the values usually 
taken in LDA+U calculationa 19 i 20 i 28 (5-8 eV). 

The Hund intra-atomic exchange integral, Jh, between 
the 3d orbitals of the cobalt atom can be evaluated from 
the excitation energies between the ground state and the 
higher spin states (S = 3/2, S = 5/2) of the Co 4+ ion. 
It comes Jh — 276 meV, to be compared with the raw 
e' g -ai g exchange integral of 0.95 eV. Hund's exchange 



values are thus strongly renormalized by the screening 
effects in Nao.35Co02 — 1.3H 2 0, leading to the low-spin 
ground-state observed in the C0O2 layers. Indeed, the 
low-spin, high-spin excitation energy can be written as 
2(e eg — e e i) — 10 Jh- With a Hund's exchange of about 
1 eV, a low spin ground state would necessitate an e g -e' g 
splitting of more than 5 eV, in total disagreement with 
all experimental and theoretical results. 



IV. DISCUSSION AND CONCLUSION 

In the present work we determined the effective on-site 
and coupling parameters for the Nao.35Co02 — 1.3H 2 
compound from ab initio quantum chemical calcula- 
tions properly treating both the strong correlation ef- 
fects within the cobalt 3d shell and the screening effects 
on the effective parameters. We determined the ligand 
field splitting as well as the on-site Coulomb repulsion 
and Hund's exchange within the Co 3d orbitals. As far 
as the interactions between two cobalt atoms are con- 
cerned, we evaluated both the effective transfer integrals 
between the t 2g orbitals as well as the effective exchange. 
It is noticeable that, the ligand field splitting between 
the aig and e' g orbitals resulting from the splitting of the 
t 2g orbitals is of the same order of magnitude (and even a 
little weaker) as the transfer integrals between two near- 
est neighbor Co atoms. In our opinion, both the a\ g and 
the e' g orbitals should be taken into account in a proper 
description of this system. 
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